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ABSTRACT 


\l 

Some  elastic  materials  are  capable  of  sustaining  finite 
equilibrium  deformations  with  discontinuous  strains.  Boundary- 
value  problems  for  such  "unstable"  elastic  materials  often 
possess  an  infinity  of  solutions,  suggesting  that  the  theory 
suffers  from  a  constitutive  deficiency.  Zn  the  setting  of  the 
one-dimensional  theory  of  bars  in  tension,  the  present  paper 
explores  the  consequences  of  supplementing  the  theory  with 
further  constitutive  information.  This  additional  information 
pertains  to  the  surface  of  strain  discontinuity  and  consists  of  a 
"kinetic  relation"  and  a  criterion  for  the  "initiation"  of  such  a 
surface.  We  show  that  the  quasi-static  response  of  the  bar  to  a 
prescribed  force  history  is  then  fully  determined.  In 
particular,  we  observe  how  unstable  elastic  meterials  can  be  used 
to  model  macroscopic  behavior  similar  to  that  associated  with 
viscopl ast ici ty . 
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1.  Introduction 

Bodies  composed  of  certain  types  of  homogeneous  elastic  materials  can 
be  finitely  deformed  to  equilibrium  states  in  which  displacement  gradients, 
strains  and  stresses  suffer  jump  discontinuities  across  special  surfaces. 

Elastostatic  fields  of  this  kind  arise,  for  example,  in  continuum  mechani-  ' 

cal  treatments  of  stress  - induced  phase  transformations  in  solids  [1,2]. 

Vhen  such  jumps  in  displacement  gradient  occur  during  quasi-static, 
isothermal  motions,  the  balance  between  the  rate  of  increase  of  stored 
energy  and  the  rate  of  work  of  external  forces  associated  with  convention-  , 

ally  smooth  deformations  of  elastic  bodies  no  longer  holds.  This  balance  is 
replaced  by  one  which  includes  an  additional  effect  that  may  be  interpreted 
as  the  rate  of  work  of  a  fictitious  "driving  traction"  acting  on  the  moving 
surface  of  discontinuity  [3].  The  driving  traction  is  formally  related  to  ' 

the  notion  of  a  "force  on  a  defect"  introduced  by  Eshelby  [4]  and  discussed 
by  Rice  [5] .  f 

The  altered  energetics  of  finite  elastostatic  fields  involving  strain 
jumps  suggest  that  such  fields  might  be  used  to  model  certain  types  of  dis¬ 
sipative  behavior  in  solids.  The  circumstances,  in  fact,  are  reminiscent  in 
some  respects  of  those  present  in  the  classical  theory  of  flows  of  ideal 
fluids  in  which  shock  waves  are  present.  In  the  latter  subject,  shocks  ] 

account  in  an  idealized  way  for  the  neglected  dissipative  effects  of  visco¬ 
sity  and  heat  conduction;  see  p.  322  of  [6].  Because  of  this  similarity,  we 
refer  to  surfaces  bearing  jump  discontinuities  in  the  displacement  gradient 
in  an  elastostatic  field  as  "equilibrium  shocks".  ] 


Not  all  elastic  materials  are  capable  of  sustaining  deformations  with 
equilibrium  shocks.  Those  that  do  have  this  capability  are  sometimes  called 
unstable  materials;  they  necessarily  lead  to  differential  equations  of 
equilibrium  that  fail  to  remain  elliptic  at  all  deformations  [7].  This  in 
turn  leads  to  a  massive  loss  of  uniqueness  of  solution  for  the  boundary 
value  problems  of  elastostatics ,  suggesting  the  need  for  additional  consti¬ 
tutive  assumptions  that  will  select  from  among  the  many  possible  equili¬ 
brium  states  that  one  which  is  preferred  by  the  body.  One  such  additional 
constitutive  postulate  asserts  that  the  material  is  conservative  at  all  of 
its  particles,  including  those  on  shocks,  so  that  the  body  prefers  that 
equilibrium  state  which  renders  the  appropriate  energy  functional  an  abso- 


! 

i 
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lute  minimum.  With  this  assumption  in  force,  the  driving  traction  acting  on 
any  equilibrium  shock  necessarily  vanishes  [8,9],  the  conventional  balance 
between  work  and  energy  is  preserved,  and  no  dissipation  takes  place.  In 
this  conservative  setting,  elastostatic  fields  with  shocks  have  recently 
received  much  analytical  attention;  see,  for  example,  [1] , [2] , [9-16] . 

In  two  recent  papers  [17,18],  we  have  discussed  an  example  in  order  to 
illustrate  an  alternative  constitutive  postulate.  The  problem  treated  in 
[17,18]  involves  a  finite,  plane  deformation  of  an  infinite  medium  contain¬ 
ing  a  circular  cavity.  A  uniform  circumferential  traction  is  applied  to  the 
cavity  wall,  and  the  displacement  is  required  to  vanish  at  infinity.  For 
the  class  of  incompressible,  isotropic  elastic  materials  considered,  the 
resulting  twisting  deformation  may  exhibit  a  circular  equilibrium  shock 
concentric  with  the  cavity.  In  quasi-static  motions  of  the  body  involving 
such  equilibrium  states,  the  relationship  between  the  applied  torque  and 
the  twist  at  the  cavity  wall  -  i.e. ,  the  macroscopic  response  -  is  in  gen¬ 
eral  hysteretic.  We  show  in  [17,18]  that  if  a  certain  maximum- dissipation 
postulate  is  used  as  the  supplementary  constitutive  assumption,  the  macros¬ 
copic  response  mimics  that  associated  with  rate - independent  elastic-plastic 
behavior . 

Our  purpose  in  the  present  paper  is  to  discuss  supplementary  constitu¬ 
tive  models  for  elastic  fields  capable  of  sustaining  equilibrium  shocks  in 
more  generality.  The  principal  new  feature  introduced  here  is  a  "kinetic 
relation"  analogous  to  those  arising  in  microstructural  models  of  plastic 
behavior  formulated  in  terms  of  internal  variables  [5,19,20].  In  our  cir¬ 
cumstances,  this  relation  takes  the  form  of  a  constitutive  law  connecting 
the  driving  traction  acting  on  a  moving  shock  with  the  shock  velocity  dur¬ 
ing  a  quasi-static,  isothermal  motion.  We  show  that  appropriate  choices  of 
the  kinetic  relation  lead  to  visco-plastic  macroscopic  response,  and  we 
recover  conservative  (minimum -energy)  response  as  well  as  rate- independent 
elastic-plastic  behavior  as  special  or  limiting  cases. 

When  there  is  an  equilibrium  shock  in  the  body,  the  kinetic  law  governs 
its  evolution.  However  a  separate  criterion  --an  initiation  or  nucleation 
criterion  --  is  required  in  order  to  signal  the  initial  appearance  of  the 
shock.  This  too  will  be  discussed  in  the  following. 


For  simplicity,  we  work  here  in  the  context  of  a  one -dimensional  model 
for  extensional  deformations  of  an  elastic  bar.  Our  setting  is  thus  essen¬ 
tially  that  of  Ericksen  [10]  in  his  discussion  of  one -dimensional  deforma¬ 
tions  with  strain-jumps ,  except  that  we  consider  bars  whose  cross-sectional 
area  varies  with  position.  The  special  case  of  the  uniform  bar  turns  out  to 
be  exceptional  in  certain  important  respects. 

After  introducing  in  the  following  section  the  class  of  elastic  materi¬ 
als  to  be  considered,  we  investigate  equilibrium  states  with  a  single  shock 
in  Section  3.  Sections  4  and  5  are  concerned  with  the  energetics  of  quasi¬ 
static,  isothermal  motions  of  the  bar  and  the  admissibility  of  such  motions 
according  to  the  second  law  of  thermodynamics.  We  introduce  the  notion  of  a 
kinetic  relation  as  well  as  a  shock  initiation  criterion  in  Section  6,  and 
in  Section  7  present  examples  to  illustrate  the  possibilities  offered  by 
the  theory. 

An  approach  of  the  type  put  forward  here  may  have  application  to  the 
modeling  of  the  mechanical  response  of  shape-memory  alloys  [21],  to  contin¬ 
uum  descriptions  of  the  effect  of  the  presence  of  a  "damaged  phase"  on  the 
behavior  of  solids,  and  to  transformation  toughening  in  ceramic  composites 
[22], 


2.  Preliminaries 

Consider  a  bar  composed  of  a  homogeneous  elastic  material,  which  in  its 
reference  configuration  occupies  the  interval  [0,L].  Let  x  denote  the  coor¬ 
dinate  of  a  generic  point  of  the  bar  in  this  configuration.  If  the  refer¬ 
ence  cross-sectional  area  of  the  bar  at  x  is  A(x)  >0,  it  is  assumed  that 
A  €  C*[0,L].  We  also  assume  that  A(x)  increases  monotonically  with  x,  so 
that 


A'(x)  >  0,  0  s  x  S  L.  (2.1) 

We  shall  show  later  that  the  special  case  of  the  uniform  bar  (A  ’(x)  -  0)  is 
exceptional  in  certain  respects;  it  is  temporarily  excluded  from  consider¬ 
ation. 


A  deformation  of  the  bar  is  characterized  by  an  invertible  mapping 
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y  -  x  +  u(x) ,  0  ^  x  5  L,  (2.2) 

which  subjects  the  particle  at  x  to  a  displacement  u  and  carries  it  to  a 
new  location  y.  There  is  no  loss  of  generality  in  taking  the  left-hand  end 
of  the  bar  to  be  fixed;  if  S  denotes  the  elongation  of  the  bar, 

u(L)  -  6,  u(0)  -  0.  (2.3) 

It  will  be  necessary  in  the  following  to  consider  displacement  fields  which 
are  less  than  classically  smooth,  and  accordingly  we  allow  for  the  possi¬ 
bility  that,  although  u  is  continuous  on  [ 0 , L] ,  there  is  a  number  s  6  [0,L] 
such  that  (i)  u  is  continuously  differentiable  on  [0,s]  +  [s,L],  (ii)  u  is 
twice  continuously  differentiable  on  (0,s)  +  (s,L),  and  (iii)  u'  suffers  a 
finite  jump  discontinuity  across  x  -  s.  The  strain  e  at  a  particle  x  s  is 
defined  by 

e(x)  -  u' (x)  >  -1,  0  <  x  <  L,  x  s;  (2.4) 

the  first  inequality  in  (2.4)  assures  the  invertibility  of  the  mapping 

(2.2). 

Let  o(x)  be  the  nominal  stress  field  in  the  bar.  Equilibrium  in  the 
absence  of  body  forces  requires 

a(x)A(x)  -  F  -  constant,  0  <  x  <  L;  (2.5) 

F  denotes  the  force  in  the  bar.  Clearly,  a  e  C^[0,L]. 

The  material  is  characterized  by  an  elastic  potential  W  whose  value 
is  the  strain  energy  per  unit  reference  volume.  We  assume  that  W  is  defined 
on  (-1,*)  and  that  it  is  twice  continuously  differentiable  there.  The 
stress  response  function  of  the  bar  <x(e)  is  given  in  terms  of  W  by 

£(«)  -  W'(€),  -1  <  £  <  (2.6) 

so  that  by  (2.5),  the  stress  at  x  is 
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S<€(x))  -  F/A(x),  o  <  x  <  L,  x  *  s.  (2.7) 

If  F  is  given,  the  force  problem  consists  of  finding  a  displacement  field 
u  of  the  requisite  smoothness  conforming  to  (2.7),  (2.4),  (2.3)2-  If  5  is 
given,  the  elongation  problem  requires  the  determination  of  a  constant  F 
and  a  displacement  field  u  satisfying  (2.7),  (2.4)  and  (2.3).  We  shall  be 
concerned  only  with  the  force  problem. 

From  a  thermodynamic  viewpoint,  the  present  analysis  assumes  that  con- 
ditons  are  isothermal.  The  elastic  potential  W  coincides  with  the  Helmholtz 
free  energy  of  the  material  at  the  given  temperature,  while  the  associated 
Gibbs  free  energy  G  expressed  in  terms  of  strain  is 


G(«)  -  W(e)  -  a(e)e  ,  -1  <  £  <  ®  .  (2.8) 


In  this  paper  we  restrict  attention  to  materials  whose  stress  response 
function  a(e)  first  increases  with  increasing  e,  then  decreases,  and 
finally  increases  again;  see  Figure  1.  Specifically  we  suppose  that  there 
are  positive  numbers  «M  and  «m  such  that  o'(eM)  -  a'(em)  -  0,  a'(e)  >  0  for 
-1  <  «  <  v'(«)  <  0  for  eM  <  «  <  em,  and  o'(e)  >  0  for  em  <  e 

Moreover , 


aM  "  a<eM>  >  °. 


£( «m)  >  0,  <7(®)  -  ®,  S(-l)  -  -®. 


(2.9) 


Note  that  these  materials  are  of  "Baker- Ericksen  type"  in  the  sense  that 
o(e)e  >  0  for  e**0 .  For  our  purposes,  it  is  sufficient  to  consider  only 
tensile  stresses,  so  we  restrict  attention  henceforth  to  a  >  0. 


Although  o(e)  is  not  invertible  on  (-1,®),  its  restrictions  to  certain 
its  of  this  interval  d 
the  analysis  to  follow. 


subsets  of  this  interval  do  have  inverses,  and  these  play  a  major  role  in 

*et  «]_,  *3  be  the  functions  inverse  to  the 
restrictions  of  o(e)  to  the  respective  intervals  (-l,eM],  [«M*eml>  and 
[ «m.®) ;  these  inverse  functions  are  defined  on  (-®,oM],  and  [am>®) 
respectively.  Each  function  «£  is  continuous  on  its  domain  of  definition, 
and  is  continuously  differentiable  on  the  interior  of  that  domain.  Finally, 
let  ac  be  the  unique  number  in  the  interval  (crm , crM)  for  which  the  two 
shaded  regions  in  Figure  1  have  equal  areas.  In  terms  of  the  Gibbs  free 
energy, 
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G(<3<*c>>  -  G(£l(ac)); 


(2.10) 


ac  is  the 


of  the  material . 


If  «(x)  is  a  solution  of  (2.7)  of  the  requisite  smoothness,  it  follows 
with  the  help  of  (2.1)  that  e(x)  *  em,£M  f°r  a H  x  in  (0,s)+(s,L),  and  that 
in  fact 


p(F/A(x) ) , 
q(F/A(x) ) , 


0  <  x  <  s , 
s  <  x  <  L, 


(3.1) 


where  p,  q  -  1,  2,  or  3.  Moreover,  for  0  <  x  <  s,  F/A(x)  must  lie  in  the 
domain  of  e'p,  while  for  s  <  x  5  L  it  must  lie  in  the  domain  of  tq.  On  using 
(2.1)  and  the  definition  of  the  inverse  functions  ,  one  can  show  that 
this  is  equivalent  to  requiring 


(  S  ,  F)  €  3r,a  , 


(3.2) 


where  the  3pq's  are  sets  in  the  (s,F) -plane  defined  as  follows: 


-  ( (s  » F) 

-  {  (s , F) 

-  {(s.F) 

-  (  (s , F) 

-  {(s.F) 

-  {(s.F) 

-  { (s , F) 

-  {  (s , F) 

-  l(s,F) 


0  <  s  <  L,  F  <  oMAm  ) , 

0  <  s  <  L,  <  F  <  aMA„,  ) , 

0  <  s  <  L,  omAM  <  F  } , 

0  <  s  <  L,  omAfl  <  F  <  ffj^Ajj,)  , 

0  <  s  <  L,  <7mA(s)  <  F  <  ), 

0  <  s  <  L,  <7mAM  <  F  <  a^}  , 

0  <  s  <  L,  CTjnA^  <  F  <  <7MA(s)  }  , 

0  <  s  <  L,  amA(s)  <  F  <  <7MA(s)  }, 

0  <  s  <  L,  <  F  <  a^Ajj,) . 


(3.3) 


Conversely,  if  (s.F)  €  3pq  for  some  p.q,  then  (3.1)  is  a  solution  of  (2.7). 
Observe  that  the  sets  322.  312>  323>  313  are  non-empty  if  and  only  if  the 
constitutive  law  and  the  taper  of  the  bar  are  such  that 


ffm^M  ^  • 


(3.4) 
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For  a  given  material,  (3.4)  will  certainly  be  valid  if  the  taper  in  the  bar 
is  slight  enough;  we  assume  throughout  that  (3.4)  holds. 

For  a  given  F,  all  solutions  of  the  force  problem  (2.7),  (2.4),  (2.3)2 

may  now  be  found  by  integrating  (3.1).  They  are 

u(x)  -  Uoa(x;F,s) ,  0  2;  x  <  L,  (s,F)  e  3_a,  p.q  -  1,2,3,  (3.5) 


where 


?p(F/A(0)  d€, 


0  <  x  <  s , 


Upq(x;F,s)  - 


(3.6) 

f  fX 

tp(F/A(0)  d£  +  tq(F/A(£))  d*.  s  <  x  <  L. 

A  s 


For  p  -  q,  (3.5),  (3.6)  yield  the  special  solutions 


u(x)  -  Up(x;F)  -  Upp(x; F , s) 


(F/A(£))  d?, 


(3.7) 


which  are  independent  of  s  and  classically  smooth.  On  the  other  hand,  for 
p  **  q,  (3.5),  (3.6)  provide  six  one-parameter  families  of  solutions  to  the 
force  problem,  with  parameter  s.  The  strain  for  0  <  x  <  s  is  associated 
with  the  pth  branch  of  the  stress -strain  curve,  while  that  for  s  <  x  <  L  is 
associated  with  the  qth  branch;  the  discontinuity  at  x  -  s  is  called  a 
" (p.q) -shock" . 

According  to  (3.3),  (3.4),  there  exists  at  least  one  solution  u(x)  to 
the  force  problem  corresponding  to  every  given  value  of  F.  If  either 
0  <  F  S  or  F  ^  o^A^,  this  solution  is  unique  and  it  is  smooth. 
However,  for  values  of  force  in  the  intermediate  range  <  F  <  o^A^, 
there  are  infinitely  many  solutions . 

Observe  from  (3.6),  (3.7)  that  as  the  shock  recedes  to  either  one  of 


jrjr 
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the  two  ends  of  the  bar,  each  weak  solution  "merges"  with  one  of  the  smooth 
solutions : 


Upq(x;F,s)  -  Uq(x;F), 
upq(x; F 1  s)  -  Up(x;F).  . 


•(3.8) 


This  suggests  that  the  definitions  of  all  of  the  Upq's,  as  functions  of  s, 
be  extended  to  s  —  0  and  s  —  L  by  setting 


Upq(x;F,0) 

-  Uq(x;F),  ' 

► 

(3.9) 

Upq(x;F,L) 

-  Up(x;F).  . 

Finally ,  one  can 

verify  that  Upq(*;F,s1) 

and  Umn(*;F,s2)  are  distinct 

whenever  (p,q)  *  (m,n) ,  0  <  s^  <  L,  and  0  <  s2  <  L: 

upq(*;F,Si)  Umn(*;F,s2)  if  (p,q)  *  (m,n),  0  <  s^  <  L,  0  <  s2  <  L, 

(Sl,F)  G  3pq,  (s2,F)  €  3mn  .  (3.10) 

We  turn  next  to  the  relation  between  the  force  F  and  the  elongation  6 , 
which  we  call  the  macroscopic  response  of  the  bar.  By  (2.3)].,  (3.5), 
(3.6),  these  quantities  are  related  by 

6-Apq(F,s),  (s,F)  G  3pq,  p,q  -  1,2,3,  (3.11) 

where 

Apq(F,s)  -  Upq(L;F,s),  (s,F)G3pq,  p.q-  1,2,3.  (3.12) 

It  can  be  shown  that,  if  p  *  q,  Apq(F,s)  is  a  monotonic  function  of  s  for 
each  fixed  F.  The  macroscopic  response  corresponding  to  any  one  of  the 
smooth  solutions  is  independent  of  s: 

S  -  App(F.s)  -  Ap(F).  (3.13) 

For  each  (p,q),  (3.11)  maps  the  set  3pq  of  the  (s,F) -plane  onto  a  set 
<Spq  in  the  (5, F) -plane: 
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-  {  (S, F) 

-  (  (S,F) 

-  {  (S, F) 

-  {  (8, F) 

-  {  (S,F) 

-  {  (5, F) 

-  (  (S,F) 

-  (  («,F) 

-  (  (8, F) 


5  -  A]_(F)  ,  F  <  }, 

S  -  A2(F) ,  <7,^  <  F  <  £7MAm  } , 

S  -  43(F),  <  F  ), 

Ai2(F,L)  <  8  <  Aj^F.O),  ^  F  <  a^Au,  )» 

A2i(F,0)  <  S  <  A2i(F,L),  ffmACs)  <  F  <  }, 

a23^f>l)  ^  8  <  A23(F,0),  <7mAfl  <  F  <  a^Am  ), 

A32(F,0)  <  6  <  A32(F,L),  ^mAM  ^  F  <  ctmA(s)  }, 

A31 (F, 0)  <  5  <  A3^(F,L),  CTmA(s)  <  F  <  aMA(s)  }, 
A13(F,L)  <  5  <  A13(F,0),  amAM  <  F  <  }. 


(3.14) 


Sketches  of  the  sets  Gpq  are  shown  in  Figure  2,  Observe  that  6^3  and  C33 
are  curves  with  positive  slope,  while  G22  a  curve  with  negative  slope. 
Note  also  that  6^3,  G22  and  C33  are  not  connected.  The  sets  Gpq,  p*q, 
correspond  to  various  regions  linking  these  curves.  The  dashed  curves  in 
these  regions  are  curves  of  constant  s.  For  p  q,  the  mappings 
3pq  -*  Gpq  are  one-to-one;  this  is  obviously  not  the  case  when  p  -  q. 

In  summary,  for  sufficiently  small  and  sufficiently  large  values  of  the 
force  F,  the  force  problem  has  a  unique  solution;  this  solution  happens  to 
be  smooth.  However,  for  intermediate  values  of  F,  we  encounter  a  major 
breakdown  of  uniqueness.  In  fact,  in  the  intermediate  range  of  F,  there  are 
multiple  solutions  even  if  the  pair  (5,F)  is  prescribed. 


4 .  Dissipation,  shock  driving  traction,  admissibili 


We  now  turn  our  attention  to  quasi-static  motions  of  the  bar  in  which, 
at  each  instant  t,  the  displacement  field  u(*,t)  corresponds  to  one  of  the 
equilibrium  states  constructed  in  the  preceding  section.  Let  F(t),  tg  <  t 
<  tj.,  be  a  given  continuous,  piecewise  continuously  differentiable  force 
history.  Suppose  first  that  F(t)  <  for  all  t  in  [tQ,t]J.  Then  by 

(3.3),  (3.7),  u(x,t)  is  necessarily  given  by  the  smooth  field 


u(x , t )  -  U1(x;F(t)) , 


0  <  x  <  L,  t0  <  t  <  t1( 


associated  with  the  first  branch  of  the  stress  -  strain  curve.  Next,  suppose 
that  F(t)  >  OflAfl  for  all  t.  Then  (3.3),  (3.7)  yield 


u(x,t)  -  U3(x;F(t)) , 


0  —  x  <  L,  t0  <  t  <  tlt 


(4.2) 


and  again  the  field  is  smooth.  Finally,  assume  that  <  F(t)  <  £7maM  f°r 

all  t  in  [tQ.t].]-  Then  u(x,t)  must  have  the  form 


u(x,t)  -  Up(t)q(t)(x;F(t) ,s(t)) ,  0  <  x  <  L,  t0  <  t  <  tx. 


To  begin  with,  we  assume  that 


(4.3) 


(i)  p(t)  and  q(t)  are  piecewise  constant  on  [tQ.tjJ,  each  taking  one  of 
the  values  1,2,3  there,  with  p(t)  #  q(t); 

(ii)  s(t)  is  piecewise  continuous  on  [tQ.t^]; 

(iii)  (s(t).F(t))  e  3p(C)q(t)  for  t0  <  t  <  t]_. 

The  requirement  p(t)  *  q(t)  does  not  preclude  the  occurrence  of  smooth 
fields  for  forces  F(t)  in  this  intermediate  range;  such  fields  occur  when 
either  s(t)  -  0  or  s(t)  -  L. 

For  quasi-static  motions  of  the  form  (4.1)  or  (4.2),  the  assumed 
smoothness  of  F(t)  and  the  representations  (3.7)  guarantee  that  u(x,*)  is 
continuous  and  piecewise  continuously  differentiable  on  [tg,t^]  for  each  x. 
In  order  to  discuss  the  more  complicated  issue  of  the  smoothness  in  time  of 
motions  described  by  (4.3),  it  is  convenient  first  to  introduce  the  notion 
of  a  transition  instant:  an  instant  t*  e  (tg,t3)  is  a  transition  instant  if 
(p(t*-) ,q(t*-) )  r*  (p(t*+) , q(t*+) ) .  At  a  transition  instant,  the  branches  of 
the  underlying  stress-strain  curve  involved  in  the  deformation  (4.3)  change 
from  the  p(t-)th  and  q(t-)th  to  p(t+)th  and  q(t+)th. 

It  is  natural  to  require  that  u(x,»)  as  given  by  (4.3)  also  be  continu¬ 
ous  on  [tg.tjJ  for  every  x,  0  <  x  <  L.  Let  t*  be  in  (tg.t^),  and  assume 
that  u(x,*)  is  continuous  at  t*.  Suppose  first  that  t*  is  not  a  transition 
instant.  Then  p(t-)  -  p(t+)  ■  p,  q(t-)  -  q(t+)  ■  q,  so  from  (4.3), 


u(x,t*+)-u(x,t*-)  -  Upq(x;F(t*) ,s(t*+))  -  Upq(x ; F( t*) , s ( t*- ) )  -  0  (4.4) 

for  every  x  in  [0,L],  Now  p  +  q,  and  the  definition  (3.6)  of  Upq  shows  that 
(4.4)  cannot  hold  under  this  circumstance  unless 


A.  «T.  m  .  A,  w*.  A  f  + .  c. ' 
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s(t*+)  -  s(t*-)  if  t*  is  not  a  transition  instant.  (4.5) 

Thus  s(t)  is  necessarily  continuous  at  all  times  except  possibly  at  transi¬ 
tion  instants.  Now  suppose  that  t*  is  a  transition  instant.  Let  p(t*-)  -  p, 
q(t*-)  ■  q,  p(t*+)  ■  m,  q(t*+)  ■  n;  by  (4.3),  continuity  of  u(x,*)  at  t-t* 
then  implies  that 

Upq(x;F(t*) ,s(t*+))  -  Umn(x;F(t*) , s(t*-))  -  0  (4.6) 

for  all  x  in  [0,L],  Since  t*  is  a  transition  instant,  ( p,q )*  (m,n),  and 
(3.10)  shows  that  (4.6)  cannot  hold  unless  at  least  one  of  the  numbers 
s(t*+) ,  s(t*-)  takes  either  the  value  0  or  the  value  L.  More  detailed 
examination  of  (4.6)  shows  that  one  of  the  following  four  mutually  exclu¬ 
sive  possibilities  must  hold: 


s(t*+)  -  s(t*-)  -  0,  and  q(t*+)  -  q(t*-),  if  (4.7a) 
s(t*+)  -  s(t*-)  -  L,  and  p(t*+)  -  p(t*-)  t*  is  a  trans-  (4.7b) 
s(t*+)  -  L,  s(t*-)  -  0,  and  p(t*+)  -  q(t*-),  ition  instant.  (4.7c) 
s(t*+)  -  0,  s(t*-)  -  L,  and  p(t*-)  -  q(t*+),J  (4.7d) 


Thus  (4.5),  (4.7)  are  necessary  for  the  continuity  of  u(x,*)  at  an 
instant  t*  in  (tg.t^);  if  either  t*  -  tg  or  t*  -  t]_,  (4.5),  (4.7)  continue 
to  be  necessary,  provided  the  appropriate  +  or  -  is  deleted  in  the  argu¬ 
ments  of  s,  p  and  q.  Moreover,  one  can  show  that,  in  the  presence  of  the 
assumed  smoothness  of  F(t) ,  (4.5),  (4.7)  (or  their  modified  versions  when 
t*  is  an  end-point  of  [tg,t^])  are  sufficient  for  the  continuity  of  u(x,*) 
at  t  ■  t*  as  well. 

The  argument  above  shows  that  discontinuities  in  s(t)  can  only  occur  at 
transition  instants  t*;  if  there  is  such  a  discontinuity,  the  shock 
x  -  s(t)  recedes  to  one  end  of  the  bar  as  t  -*  t*-  and  then  advances  into 
the  bar  from  the  other  end  as  t  increases  from  t*.  When  s(t)  is  continuous 
at  a  transition  instant  t*,  necessarily  s(t*)  -  0  or  s(t*)  -  L.  Thus  a 
transition  from  a  discontinuous  strain  field  involving  branches  p  and  q  of 
the  stress  -  strain  curve  to  one  involving  branches  m  and  n,  with  (p,q)  ** 
(m,n),  always  takes  place  through  a  smooth  field.  A  further  consequence  of 
the  conditions  (4.5),  (4.7)  is  that  a  shock  cannot  emerge  instantaneously 
from  a  smooth  field  at  an  interior  point  of  the  bar.  Observe  that  these 


restrictions  on  the  motion  x  —  s(t)  of  the  shock  arise  from  purely  kine¬ 
matic  requirements,  together  with  the  assumption  that  the  bar  is  strictly 
monotonically  tapered.  Further  restrictions  on  the  shock  motion  will  arise 
later. 

Finally,  we  require  that  s(t)  should  be  piecewise  continuously  differ¬ 
entiable  between  every  pair  of  successive  transition  instants.  This  will 
assure  that  u(x,*)  is  piecewise  continuously  differentiable  on  [tQ.tjJ. 
A  regular  instant  during  a  quasi-static  motion  is  a  time  t  at  which  F(t) 
exists  and,  if  the  motion  is  of  the  form  (4.3),  s(t)  exists  and  p(t)  and 
q(t)  are  continuous. 

Figure  3  describes  an  example  in  the  (x.t) -plane  for  which  the  shock 
history  involves  transition  instants  of  each  of  the  four  kinds  listed  in 
(4.7a-d).  The  encircled  numbers  in  Figure  3  refer  to  the  branches  of  the 
stress-strain  curve  appropriate  to  the  two  sides  of  the  shock  at  various 
times . 

The  elementary  quasi-static  motions  (4.1), (4. 2)  and  (4.3)  may  be  linked 
together  on  successive  time  intervals  [  t0,  t]J  ,  [  tj_,  t2]  ,  and  [t2,t3]  to  form 
a  compound  quasi-static  motion  on  [t0,t3],  provided  the  resulting  displa¬ 
cement  u(x , • )  is  continuous  on  [tg,t3]  for  every  x. 

Next  we  consider  the  energetics  of  a  quasi-static  motion.  The  total 
strain  energy  stored  in  the  bar  in  an  equilibrium  state  with  displacement 
Upq(x; F , s)  is 


E 


pq 


CF.S) 


|  W(fp(F/A(x)))  A(x)  dx  +  J  W(eq(F/A(x)))  A(x)  dx, 

P-q  -  1 '2-3.  (s.F)  6  3pq.  (4.8) 


During  a  quasi-static  motion  of  the  form  (4.3),  the  energy  stored  in  the 
bar  at  time  t  is 


E(t)  -  Ep(t)q(t)(F(t)  ,s(t)) ,  t0  <  t  <  tj_. 


(4.9) 


At  a  regular  instant  during  this  motion,  we  define  the  rate  of  dissiDation 


-.U-- 


d(t)  to  be  the  difference  between  the  rate  of  external  work  and  the  rate  of 
increase  of  stored  energy: 


d(t)  -  F(t)5(t)  -  E(t), 


(4.10) 


where,  by  (3.11),  the  elongation  5(t)  of  the  bar  is 


S(t)  -  Ap(t)q(t)(F(t),s(t)),  tQ  <  t  <  tx, 


(4.11) 


and  Apq  is  given  by  (3.12).  A  direct  calculation  based  on  (4. 8) -(4. 11) 
yields 


d(t)  -  [ [W(e)  -  o(«)f] ]+A(s(t))s(t) , 


(4.12) 


for  each  regular  instant  during  the  motion.  Here  [[g]]  denotes  the  jump  at 
x  -  s(t)  of  the  function  g(x) :  [ [ g ] J  -  g(s(t)+)  -  g(s(t)-).  It  follows  from 
(4.12)  that  d(t)  -  0  if  the  motion  happens  to  be  smooth  at  time  t,  so  that 
all  jumps  in  (4.12)  vanish.  In  general,  however,  d(t)  *  0  for  motions  of 
the  type  (4.3).  For  motions  of  the  type  (4.1)  or  (4.2),  the  dissipation 
rate  of  course  vanishes  identically. 


Let  f(t)  be  defined  by 


f(t)  -  [  (W(  e )  -  o(0«3].  .  tQ  <  t  <  tx, 


(4.13) 


and  by  f(t)  -  0  for  motions  of  the  type  (4.1)  or  (4. 2). Since  by  (4.10), 

(4.12),  (4.13),  at  a  regular  instant 


E  -  FS  +  (-f  A(s) )  s. 


(4.14) 


it  follows  that  one  can  view  -f(t)  as  a  traction  applied  by  the  shock  on 
the  bar,  or  equivalently,  +f(t)  as  a  traction  applied  by  the  bar  on  the 
shock.  The  general  notion  of  a  "thermodynamic  force"  or  "driving  force”  on 
a  "defect"  was  introduced  by  Eshelby  [4]  .  Equation  (4.13)  is  a  special 
case  of  a  formula  for  the  force  on  a  phase  boundary  given  in  [4]  and  dis¬ 
cussed  by  Rice  [5];  see  also  [3].  We  refer  to  f  as  the  "shock  driving  trac- 
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tion" .  Observe  from  (2.8),  (4.12),  that  the  shock  driving  traction  coin¬ 
cides  with  the  jump  in  the  Gibbs  free  energy  across  the  shock: 


f(t)  -  [[G(«(x,t))]]_. 


(4.15) 


On  using  (2.10),  one  obtains  from  (4.15)  the  following  expression  for  the 
shock  driving  traction: 


f(t)  -  fp(t)q(t)(ff(t)) 


(4.16) 


where  <T(t)  is  the  nominal  stress  at  the  shock, 


a(t)  -  F(t)/A(s(t)), 


(4.17) 


and  the  functions  fD_  are  determined  by  the  material;  they  are  given  by 


ff(e)  d« 


-  <7{«q(a)  -  ep(o)},  C7m  <  a  <  <7M,  p,q  -  1,2,3.  (4.18) 


The  following  properties  of  the  functions  fpq  will  prove  to  be  useful. 
First,  note  from  (4.18)  that 


fpq  (  o )  "  *  ^qp  ^  a  ^  ’ 


om  <  a  <  p,q  -  1,2,3. 


(4.19) 


Second,  differentiating  (4.18)  yields 


>  0,  am  <  a  <  oM,  p  >  q. 


f'VQ(o)  “  ep(«7)  '  eq(ff)  1  “  °-  am  ^  a  ^  CTM>  P  “  9-  (4-20) 


<  0,  am  <  a  <  aM,  p  <  q. 


Next,  on  recalling  the  properties  of  the  inverse  functions  e^,  one  can 
readily  verify  from  (4,18)  that 


f31<ffm)  <  °.  f3l(CTc)  “  °>  f31<aM> 

f2l(ffM)  “  °>  ^32 ^m)  ” 


>0.  J 


(4.21) 


l 
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Integration  of  (4.20)  with  the  help  of  (4.21)  then  gives  the  following  use¬ 
ful  alternative  formulas  for  fpqi 


f32(o) 


^31  (<T) 


f2l<«0 


-r 

an 

-f 


{e3('>  *  <2<r>)  dr, 


{f3(r)  -  ?x(r))  dr, 


M> 


£Tm^<7— ctM  ’ 


•  ( 


{?2(r)  -  «i(r))  dr, 


(4.22) 


the  other  fpq's  are  related  to  these  through  (4.19).  Observe  from  (4.22) 
that  f3i,  f32  and  f2^  are  monotonically  increasing  functions  on  (CTm,£7M)- 
Moreover,  f32  is  non-negative  on  its  domain  of  definition  while  f2^  is 
non-positive;  on  the  other  hand  f3*L  changes  sign  once: 

f32(cr) 

f3l(ff)  1  “  °.  o~ac,  (4.23) 

l  <  0.  o^aKa^., 

f  —  0 .  a—Ou . 

f 21 (a) 


A  quasi-static  motion  is  said  to  be  admissible  if  the  rate  of  dissipa¬ 
tion  is  non-negative: 


> 

0, 

am<a—aH > 

- 

0, 

a“am’ 

> 

0, 

ac<0^aH < 

- 

0. 

o~oc, 

< 

0, 

av&a<ac  • 

- 

0, 

o-oM, 

< 

0, 

crm^c7<<7M  • 

f(t)s(t+)  >  0, 


CQ  <  t  <  tL. 


(4.24) 


On  adapting  the  thermodynamic  arguments  given  in  [3]  to  the  present  one- 
dimensional  context,  one  can  show  that  (4.24)  is  equivalent  to  the  Clau- 
sius-Duhem  inequality  when  the  temperature  is  constant  in  space  and  time. 
Observe  that  (4.24)  holds  with  equality  if  the  shock  is  stationary  or  if 
the  shock  driving  traction  vanishes;  the  latter  alternative  occurs  if  and 
only  if  either  (i)  the  field  is  smooth  or  (ii)  (p.q)  -  (3,1)  or  (1,3)  and 


i 

) 
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the  stress  at  the  shock  is  the  Maxwell  stress.  All  motions  of  the  type 
(4.1)  or  (4.3)  are  trivially  admissible.  In  general,  (4.24)  is  to  be  viewed 
as  a  restriction  on  allowable  quasi-static  motions.  Note  from  (4.11), 
(4.24)  that  a  quasi-static  motion  is  admissible  at  an  instant  t  if  and  only 
if  the  Gibbs  free  energy  of  a  particle  at  the  shock  front  does  not  increase 
as  the  particle  crosses  the  shock.  In  the  following  section,  we  investigate 
the  consequences  of  admissibility. 


5 .  Consequences  of  admissibility. 

Consider  an  admissible  quasi-static  motion  of  the  form  (4.3)  in  which 
p(t)  ■  p  -  constant,  q(t)  ■  q  -  constant  for  tg  <  t  <  t3.  If,  for  example, 
p-3,  q— 1,  then  by  (4.20),  (4.12)  and  (4.19),  admissibility  requires  that 


s(t+)  >  0 
s(t+)  <  0 
<T(t+)s(t+)  >  0 


if  ff(t)  >  ac, 
if  o(t)  <  ac, 
if  <F( t)  -  <rc, 


(5.1) 


where  cr(t)  •  F(t)/A(s(t))  is  the  stress  at  the  shock.  The  curve  in  83^ 


given  by  F  -  ac A(s) ,  S  -  A31(acA(s) ,s)  is  called  the  3^1 -Maxwell  curve: 
points  on  this  curve  correspond  to  equilibrium  fields  in  which  the  stress 
at  the  shock  coincides  with  the  Maxwell  stress. 


Every  quasi-static  motion  determines  a  path  in  the  (5, F) -plane.  For  the 
motion  considered  above,  this  path  lies  in  the  set  G33  (see  Figure  2),  and 
-  through  (5.1)  -  admissibility  restricts  the  possible  directions  of  the 
path.  Figure  4(a)  illustrates  the  permissible  directions  of  departure  of 
such  a  path  from  various  points  in  <S33.  The  dashed  curves  in  the  figure 
represent  lines  s(t)  -  constant;  for  motions  whose  paths  lie  along  a  curve 
s  -  constant,  there  is  no  dissipation.  The  same  is  true  for  motions  whose 
paths  lie  on  the  Maxwell  curve. 


Similar  considerations  apply  to  motions  of  the  type  (4.3)  for  other 
values  of  the  pair  (p,q).  For  p-1,  q-3,  there  is  also  a  Maxwell  curve,  but 
for  the  remaining  possible  choices  of  (p,q),  this  is  not  the  case.  Figures 
4(b) -4(f)  show  the  admissible  directions  for  the  remaining  choices  of 
(p,q).  In  Figures  4(a)  and  4(b)  and  hereafter,  we  assume  that 


■VSfii 


am/ac  <  V/Am  •  ac/aM  <  Atn/^M  • 


(5.2) 


I 

i 

These  inequalities,  which  imply  (3.4),  certainly  hold  for  a  given  material 

if  the  taper  of  the  bar  is  slight  enough.  , 

i 

i 

i 

Admissible  directions  for  compound  motions  can  also  be  read  off  from  ' 

Figure  4. 

From  Figure  4(a),  it  is  clear  that  transitions  of  the  form  ^^^-*63^, 

^Sr^ll’  ®33"*®31'  ®3r*®33  are  possible;  similar  inferences  may  be  drawn 

from  Figure  4(b).  The  situation  is  different,  however,  for  admissible  ; 

motions  that  involve  branch  2  of  the  stress -strain  curve.  For  example,  Fig-  ] 

i 

ure  4(c.)  shows  that,  while  the  transition  ®21"*®11  admissible,  the  > 

reverse  transition  is  not;  similarly,  ®22"*®21  *-s  possible,  but  G2r*622  is 
not.  Parallel  conclusions  follow  from  Figures  4(d)-4(f).  In  general,  one 
can  readily  show  that  admissible  quasi -static  motions  proceed  in  such  a  way 

i 

as  to  diminish  -  or  at  least  not  increase  -  the  length  of  the  portion  of  | 

the  bar  bearing  strains  associated  with  branch  2  of  the  stress -strain 

curve .  | 

The  above  considerations  suggest  that  the  totality  of  all  equilbrium 
displacement  fields  be  divided  into  two  disjoint  parts  A  and  U  as  follows: 
let  I 


A  -  {Upq(.;F,s)  |  (s.F)  e  3pq,  p  -  1  or  3 ,  q  -  1  or  3), 

U  -  (Upq(  *  I F,  s)  |  (s  ,  F)  €  3pq ,  either  p  -  2  or  q  -  2).  ' 

Thus  the  fields  in  A  correspond  to  those  equilibrium  states  that  involve  no 
strains  associated  with  branch  two  of  the  stress  -  strain  curve  in  Figure  1; 
U  contains  all  remaining  displacement  fields.  Consider  a  compound  admis¬ 
sible  motion  whose  displacement  field  at  the  initial  instant  belongs  to  A. 
Figure  4  suggests  that  displacement  fields  for  this  motion  at  all  later 
times  must  also  belong  to  A.  Indeed,  if  at  some  later  instant  the  corre¬ 
sponding  displacement  field  were  in  U,  the  length  of  that  portion  of  the 

bar  carrying  branch-2  strains  would  necessarily  have  increased  at  some  ear¬ 

lier  time,  contradicting  the  assumed  admissibility  of  the  motion.  Thus 
states  in  U  are  not  accessible  from  states  in  A  during  an  admissible  quasi - 
static  motion;  in  particular,  a  motion  that  commences  at  the  unloaded, 
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unextended  state  F  -  5  -  0  can  never  enter  the  collection  U  of  states 
involving  strains  on  branch  2  of  the  stress -strain  curve. 

Let  d(t)  be  the  rate  of  dissipation  in  a  quasi-static  motion  at  each 
regular  instant  t;  the  total  dissipation  D  associated  with  the  motion  is 

rCl 

D  -  d(t)  dt  .  (5.4) 

*-0 

It  is  possible  to  show  that  all  admissible  quasi-static  motions 
whose  total  dissipation  is  sufficiently  large  must  ultimately  enter  the 
collection  A,  where  -  in  view  of  the  discussion  above  -  they  will 
subsequently  remain.  We  shall  not  prove  this  result  here;  see  [17,18]  for 
proofs  of  closely  related  propositions. 

One  may  thus  regard  the  collection  A  of  equilibrium  fields  as  an 
attractor  for  admissible  quasi-static  motions;  the  set  U  may  be  thought  of 
as  unobservable.  From  here  on,  we  shall  be  concerned  only  with  motions 
through  equilibrium  states  that  can  be  reached  admissibly  from  the  state  F 
-6-0;  as  a  result,  we  need  no  longer  consider  states  that  include 
strains  associated  with  branch  2  of  the  stress -strain  curve. 


In  reaching  the  conclusions  described  above,  the  fact  that  the  cross- 
section  of  the  bar  is  tapered,  rather  than  uniform,  is  crucial.  For  a  uni¬ 
form  bar,  the  implications  of  admissibility  are  much  weaker  than  those 
described  above.  The  distinction  can  be  appreciated  with  the  help  of  Figure 
5,  which  shows  the  sets  Gpq  for  the  uniform  bar.  Observe  from  the  figure 
that  the  sets  5^,  622  and  G33  corresponding  to  smooth  strain  fields  are 
now  connected,  each  to  the  next,  in  contrast  to  the  corresponding  sets  for 
the  tapered  bar  as  shown  in  Figure  2  or  Figure  4.  Thus  for  a  tapered  bar, 
one  cannot  move  quasi-statically  from  states  in  G^i  to  states  in  622  by 
utilizing  smooth  fields  alone;  such  a  transition  would  demand  the  presence 
of  fields  with  equilibrium  shocks,  which  -  by  the  discussion  above  -  is 
forbidden  when  admissibility  is  imposed  as  a  requirement.  On  the  other  hand 
Figure  5  shows  that,  for  the  uniform  bar,  one  can  construct  quasi-static 
motions  involving  only  smooth  -  and  therefore  trivially  admissible  -  fields 
that  pass  from  states  in  G^  to  states  in  G22-  Thus  while  admissibility 


£3 


forces  the  collection  U  of  equilibrium  states  involving  branch  2  to  be 
inaccessible  and  transient  in  the  presence  of  sufficient  dissipation  in  the 
case  of  the  tapered  bar,  it  does  not  deliver  the  corresponding  results  for 
the  uniform  bar. 

6 .  Kinetic  relations  and  shock  initiation. 

The  specification  of  either  the  force  history  F(t)  or  the  elongation 
history  5(t)  during  an  admissible  quasi-static  motion  is  not  sufficient  to 
determine  the  motion  uniquely  unless  the  field  is  smooth  at  each  instant. 
This  suggests  that  the  constitutive  characterization  of  the  material  must 
be  supplemented  if  the  response  is  to  be  determinate  when  equilibrium 
shocks  are  present. 

In  the  macroscopic  F-5  relations  (3 . 11) - (3 . 13)  pertaining  to  equili¬ 
brium  states  with  shocks,  the  shock  location  s  may  be  viewed  as  an  "inter¬ 
nal  variable".  Indeed,  the  formalism  used  in  internal  variable  theories  of 
plasticity  such  as  those  developed  in  [5],  [19]  and  [20]  has  a  counterpart 
in  the  present  context.  Because  we  do  not  need  this  formalism  here,  we  do 
not  discuss  it  in  detail;  it  has  been  described  in  a  related  setting  in 
[17]  and  [18].  A  common  ingredient  of  such  microstructural  theories  of 
inelastic  behavior  is  an  evolution  law,  or  "kinetic  relation",  relating  the 
time  rate  of  change  of  the  internal  variable  to  the  corresponding  "thermo¬ 
dynamic  force".  We  adopt  this  point  of  view  here  by  postulating  a  relation 
between  the  shock  driving  traction  f(t)  of  (4.13)  and  the  shock  velocity 
s(t)  that  must  hold  during  a  quasi-static  motion.  This  relation  is  deter¬ 
mined  by  the  material  and  is  regarded  as  given. 

Suppose  we  have  an  admissible  motion  of  the  form  (4.3).  Recall  the 
relation  (4.12)  between  the  shock  driving  traction  f(t)  during  such  a 
motion  and  the  stress  at  the  shock  “(t),  and  let 

Mpq  -  max  fpq(ff),  nL,q  -  min  fpq(o),  (6.1) 

be  the  maximum  and  minimum  values  of  the  material  functions  fpq  introduced 
in  (4.18).  For  each  p,q  -  1,2,3  with  p  *  q,  we  assume  there  is  a  given 
function  Vpq  defined  on  (rapq.Mpq)  x[0,L]  and  such  that,  between  successive 
transition  instants  during  the  motion,  the  kinetic  relation 
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s(t)  “  vp(t)q(t)(f<t) ,s(t)) 


(6.2) 


holds.  If  the  functions  Vpq(f,s)  are  independent  of  s,  we  say  that  the 
kinecic  relation  is  homogeneous :  we  assume  this  to  be  the  case  henceforth. 


We  now  impose  on  each  kinetic  response  function  Vpq  the  requirement 


Vpq(f)f  >  0  ,  mpq 


<  f  < 


Mpq,  0  <  s  <  L; 


(6.3) 


by  (4.24),  (6.2),  this  assures  that  all  quasi-static  motions  compatible 


with  the  kinetic  relation  are  admissible. 


We  also  require  that,  for  each  p,q  with  p  +  q, 


Vqp(f)  -  -  Vpq(-f)  for  -  Mpq 


<  f  <  - 


(6.4) 


(Note  that,  by  (4.9),  (6.1),  mqp  -  -Mpq.)  The  motivation  for  (6.4)  lies  in 
the  fact  that  the  kinetic  response  functions  Vpq  are  to  depend  only  on  the 
material  and  not  on  the  geometry  of  the  bar  under  consideration.  In  partic¬ 
ular,  they  must  apply  in  the  case  of  uniform  bars,  for  which  (3,1) -shocks 
and  (1,3) -shocks  are  mirror- images  of  one  another  if  the  stress  at  the 


shock  is  the  same  for  both. 


Since  the  shock  driving  traction  f(t)  during  the  motion  is  related  to 
the  stress  at  the  shock  eT(t)  by  (4.16),  the  kinetic  relation  relation  (6.2) 
may  be  expressed  in  terms  of  o'  instead  of  f: 


S(t)  -  vp(t)q(t)(*(t))  -  vp(t)q(t)(F(t)/A(s(t))) , 


(6.5) 


where  the  material  function  v_q  is  defined  by 


VPQ  (<7)  “  vpq(^pq(<7) )  >  <  a  <  ^M- 


(6.6) 


By  (4.16),  (4.19),  properties  (6.3),  (6.4)  of  the  Vpq ' s  yield  corresponding 
properties  of  the  Vpq's: 


i 

»« 

* 


i 
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(6.7) 

(6.8) 


vpq(CT)fpq(a)  >  0  ,  am<  a  <  an, 

Vqp(tf)  -  -  vpq(ff) ,  am  <  a  <  ctm. 

Between  two  successive  transition  instants,  p(t)  and  q(t)  are  both  con¬ 
stant:  p(t)  ■  p,  q(t)  ■  q  ;  for  a  force-controlled  motion  in  which  F(t)  is 
given,  (6.5)  is  then  a  differential  equation  governing  the  location  x  - 
s(t)  of  the  associated  (p,q)-shock. 

For  definiteness,  consider  a  program  of  loading  in  which  the  given 
force  history  F(t),  0  <  t  <  T,  begins  at  F(0)  -  0,  so  that  initially  6(0)  - 
0  as  well,  and  suppose  that  F(t)  increases  with  t  to  a  final  value  F(T)  > 
To  describe  the  possible  quasi-static  motions  associated  with  F(t), 
we  shall  trace  the  corresponding  force-elongation  histories  in  Figures 
6(a,b),  which  contain  the  information  in  Figures  4(a)  and  4(b)  pertaining 
to  equilibrium  states  with  (3,1)-  and  (1 , 3) -shocks .  (Recall  from  the  pre¬ 
ceding  section  that  states  with  either  p  -  2  or  q  -  2 ,  or  both,  can  never 
be  reached  admissibly  during  a  force  history  of  this  kind,  so  that  Figures 
4(c) -4(f)  need  not  be  consulted.)  Reference  to  Figure  6  shows  that,  as  F(t) 
increases  from  zero  to  the  value  at,  say,  time  t^ ,  the  associated 

quasi-static  motion  is  necessarily  smooth  and  of  the  form  (4.1).  During 
this  period,  the  force-elongation  relation  corresponds  to  the  curve  OA  in 
either  Figure  6(a)  or  6(b),  and  according  to  (3 . 11) - (3 . 13) ,  it  is  given  by 

«(t)  -  Ax(F(t)).  (6.9) 


As  F(t)  increases  beyond  equilibrium  states  involving  (3,1) -shocks 

become  available  (Figure  6(a)),  but  as  indicated  by  the  arrows  in  the  fig¬ 
ure,  none  can  be  reached  admissibly  until  the  time,  say  t-2,  at  which  F(t2) 
-  WqAjj.  Thus  on  the  time  interval  [  t  ,  t£  ]  .  the  motion  is  of  the  form  (4.3), 
with  u(x,t)  -  U3^(x;F(t) ,0)  -  U^3(x;F(t) ,L)  -  Ui(x;F(t)),  the  displacement 
field  remains  smooth  at  each  instant,  (6.9)  remains  in  force,  and  the  for¬ 
ce-elongation  curve  continues  along  the  arc  ABC.  For  t  >  t2  and  F(t)  > 
the  situation  changes.  Let  F(t3)  -  <7CA^;  when  t2  <  t  <  t3>  there  are 
two  alternative  possibilities,  each  consistent  with  admissibility:  either 
the  deformation  continues  to  be  smooth  at  each  instant,  with  u(x,t)  - 
U^(x;F(t)),  or  a  (3,l)-shock  is  initiated  at  the  end  x  -  0  of  the  bar  at  a 
certain  instant  t*  >  t2  and  begins  to  advance  into  the  interior  in 
accordance  with  the  kinetic  relation  (6.2)  (or  (6.5)).  In  the  former  event. 


I'/iYt'.yjy 
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(6.9)  continues  to  hold,  and  the  force  elongation  curve  is  ABCD  (Figures 
6(a,b)).  However,  if  a  (3,l)-shock  forms,  then  the  F-6  curve  departs  from 
ABCD  at  a  point  0*  (Figure  6(a)),  with 

S  -  A31(F(t),s(t)).  (6.10) 

Here  s(t)  is  the  location  of  the  shock  at  time  t  >  t*;  it  is  given  by  the 
solution  of  the  differential  equation  (6.5)  with  p-3,  q-1,  subject  to  the 
initial  condition  s(t*)  -  0.  If  a  shock  does  not  form  and  the  first  of  the 
two  above  alternatives  occurs,  the  force,  having  generated  only  smooth 
deformations,  eventually  attains  the  value  F(tj)  -  ocA^,  after  which  any 
one  of  three  mutually  exclusive  types  of  admissible  response  histories  may 
occur.  First,  the  fields  may  remain  smooth,  with  6(t)  -  A^(F(t));  second,  a 

(3.1)  -shock  may  emerge  at  x  -  0,  after  which  (6.10)  will  take  over;  third, 
it  is  now  possible  to  create  a  (1,3) -shock  at  x  -  L.  If  this  third  possi¬ 
bility  occurs  ,  say  at  time  t**,  the  force-elongation  curve  will  depart 
from  ABODE  at  0**  (Figure  6(b)),  and  the  F-6  relation  will  be 

«(t)  -  A13(F(t),s(t)),  (6.11) 

instead  of  (6.9). 

The  kinetic  relation  alone  does  not  determine  whether  or  when  a  shock 
forms  during  the  above  loading  program  or,  if  so,  whether  it  is  a 

(3.1) -shock  at  x  -  0  or  a  (l,3)-shock  at  x  •  L.  It  is  therefore  necessary 
to  establish  in  addition  a  criterion  for  shock  initiation.  If  a  (3,1) -shock 
forms  at  x  -  0  at  time  t*  (corresponding  to  the  point  0*  in  Figure  6(a)), 
then  the  strain  at  the  particle  x  -  0  will  jump  from  a  value  associated 
with  branch  1  of  the  stress-strain  curve  to  a  value  on  branch  3;  for 
brevity,  we  say  that  the  particle  has  undergone  a  transformation  from 
"phase  1"  to  "phase  3"  at  time  t*.  The  same  can  be  said  of  the  particle  x  - 
L  if  a  (1,3) -shock  emerges  from  x  -  L  at  time  t**  (point  0**  in  Figure 
6(b)).  We  now  adopt  a  specific  criterion  for  such  shock- initiating  -  or 
"spontaneous"  -  phase  1  -»  phase  3  transformations:  the  particle  x  -  x*  will 
spontaneously  transform  from  phase  1  to  phase  3  at  time  t*  if  the  stress 

<r(x*,t*)  >  and  o(«,t*)  has  a  local  maximum  at  x*.  (6.12) 
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Here  the  "transformation  stress"  is  determined  by  the  material  and 
presumed  to  be  known.  (The  alternative  to  a  spontaneous  transformation 
occurs  when  a  particle  changes  phase  due  to  the  passage  of  a  pre-existing 
shock  through  that  particle.)  For  the  reverse  transformation,  we  similarly 
postulate  that  the  particle  x*  will  spontaneously  transform  from  phase  3  to 
phase  1  at  time  t*  if 

<r(x*,t*)  <  Sjjj  and  a(*,t*)  has  a  local  minimum  at  x*,  (6.13) 

where  the  reverse  transformation  stress  5^  is  also  given.  Admissibility 
requires  that  SM  and  Eg,  satisfy 

am  <  ^m  s  ac  -  SM  <  aM  •  (6.14) 

Since  the  bar  is  monotonically  tapered,  the  maximum  stress  at  each  instant 
occurs  at  the  small  end  x  -  0,  the  minimum  at  x  -  L.  Thus  shocks  corre¬ 
sponding  to  phase  1  -♦  phase  3  transformations  can  only  be  initiated  at  x  - 
0,  and  those  corresponding  to  3  -*  1  transformations  only  at  x  -  L. 

We  remark  that  the  shock  initiation  criterion  introduced  above  can  be 
alternatively  formulated  in  terms  of  critical  values  of  shock  driving  trac¬ 
tion.  rather  than  in  terms  of  the  critical  values  SM  and  Ejjj  of  stress  at 
the  shock.  In  the  present  one -dimensional  context,  no  advantage  is  gained 
by  using  this  formulation,  so  we  omit  it.  In  higher  dimensional  settings, 
however,  it  is  likely  that  the  alternative  formulation  is  essential. 

When  applied  to  the  loading  program  described  in  detail  above,  our 
criterion  singles  out  a  definite  response  history  from  among  the  possibili¬ 
ties  listed  there:  as  t  increases  from  zero,  the  equilibrium  fields  remain 
smooth  and  the  force  elongation  relation  remains  given  by  (6.9)  until  the 
instant  t*  at  which  o(0,t*)  -  F(t*)/A1Il  -  At  time  t*,  a  (3,l)-shock 
forms  at  x  -  0.  The  evolution  of  the  shock  location  s(t)  is  then  con¬ 
trolled  by  the  differential  equation  (6.5)2  p(t)  -  3,  q(t)  -  1,  sub¬ 
ject  to  the  initial  condition  s(t*)  -  0.  Under  suitable  restrictions  on  the 
kinetic  response  function  V31  to  be  specified  in  the  following  section,  the 
associated  force-elongation  response  curve,  described  now  by  (6.10),  will 
remain  in  the  set  333  of  (3 , 1) -equilibrium  states;  it  is  the  curve  0*P* 
shown  schematically  in  Figure  6(a). 
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As  the  force  F(t)  increases  above  the  level  associated  with  the  point 
P*,  the  subsequent  response  necessarily  is  smooth  and  corresponds  to 
branch- 3  strains;  thus  the  F-6  relation  now  becomes 

S(t)  -  U3(F(t),s(t))  (6.15) 
during  the  remainder  of  the  loading  process. 

If  the  force  F(t)  is  now  decreased  monotonically  to  zero  from  its 
greatest  value  F(T) ,  the  ensuing  deformation  will  be  smooth,  and  the  force  - 
elongation  relation  will  be  given  by  (6.15)  until  the  minimum  stress  in  the 
bar  o(L,t)  reaches  the  lower  transformation  stress  2^.  At  this  instant,  the 
particle  x  -  L  undergoes  a  3  -»  1  phase  transformation,  and  a  (3,1) -shock  is 
created  at  x  -  L,  moving  leftward  into  the  bar  in  accordance  with  (6.5). 
The  F-6  response  curve,  though  again  described  by  (6.10),  will  differ  from 
its  counterpart  during  loading.  Once  this  curve  rejoins  OABC,  the  response 
during  the  remainder  of  the  unloading  process  is  that  associated  with 
branch-1  smooth  fields,  continuing  down  OABC  to  the  origin. 

Finally,  it  should  be  noted  that  kinetic  relations  whose  structure  is 
substantially  more  general  than  (6.2)  could  be  considered. 


7 .  A  special  class  of  kinetic  relations. 

We  now  consider  a  special  class  of  kinetic  response  functions  V  in 
(6.2)  (or  v  in  (6.5)).  After  investigating  some  of  their  properties,  we 
illustrate  in  detail  the  visco-plastic  nature  of  the  macroscopic  response 
of  the  bar  to  which  they  give  rise.  We  discuss  rate  -  independent  dissipative 
response  and  purely  conservative,  dissipation- free  behavior  in  the  present 
context  as  well. 


Since  we  shall  be  concerned  only  with  equilibrium  states  accessible 
through  admissible  motions  from  the  unloaded,  unextended  state  F  -  5  -  0 
we  will  not  encounter  strains  on  branch  2  of  the  stress-strain  curve.  More¬ 
over,  for  the  loading  programs  to  be  considered,  the  shock  initiation 
criterion  of  the  preceding  section  will  always  rule  out  (1 , 3) -shocks .  Thus 
we  shall  be  concerned  with  the  kinetic  relation  (6.2)  only  when  p(t)  -  3, 
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q(t)  -  1.  As  a  result,  we  shall  write  V31  -  V,  Vp^  -  v,  M31  *  M,  -  m 
throughout  the  present  section  for  convenience.  Recall  from  (4.22),  (4.23) 
and  (6.1)  that  m  <  0,  M  >  0. 


Guided  in  part  by  the  form  sometimes  ascribed  to  the  counterpart  of  the 
function  V  in  microstructural  theories  of  plasticity  [24],  we  now  make 
three  assumptions  concerning  the  form  of  V:  (i)  We  assume  that  there  are 
numbers  m*  and  M*  such  that 


m<m*<0<M*<M 


(7.1) 


V(f)  <  0  for  m  <  f  <  m*, 


V(f)  -  0  for  m*  <  f  <  M*, 


V(f)  >  0  for  M*  <  f  <  M, 


(7.2) 


noting  that  (7.2)  is  consistent  with  the  requirement  (6.3)  imposed  by 
admissibility.  (ii)  We  assume  that  V(f)  is  continuous  on  (m,M)  and  conti¬ 
nuously  differentiable  on  (m,m*]  +  [M*,M) ,  and  that  V' (f)  >  0  for  m  <  f  < 
m*  and  for  M*  <  f  <  M.  (iii)  Finally,  we  assume  that  V(f)  is  unbounded  as  f 
-♦  nH-  and  as  f  -*  M- ;  more  precisely,  we  require  that 


V(f)  <  Cm(f  -  m)'1  for  f  sufficiently  near 
V(f)  >  C[j(M  -  f)’^  for  f  sufficiently  near 


m,  1 

M,  J 


(7.3) 


for  suitable  constants  Cm  <  0  and  >  0. 


According  to  (i),  a  (3,l)-shock  will  move  in  the  +x-direction  only  if 
the  shock  driving  traction  exceeds  M*,  and  in  the  reverse  direction  only  if 
f  is  less  than  m*.  Permitting  V(f)  to  vanish  over  the  interval  [m*,M*]  will 
be  seen  later  to  allow  for  the  possibility  of  dissipation-free  unloading 
and  re-loading.  Assumption  (ii)  assures  that,  roughly  speaking,  larger 
shock  tractions  promote  greater  shock  speeds.  The  role  of  property  (7.3) 
will  become  clear  as  we  proceed.  A  schematic  sketch  of  the  graph  of  V(f)  is 
shown  in  Figure  7(a);  it  is  reminiscent  of  the  relationship  sometimes  pro¬ 
posed  between  the  driving  force  on  a  dislocation  and  dislocation  velocity 


By  (6.6),  property  (7.2)  of  V  translates  into  a  corresponding  property 


of  v: 


v(o)  <  0  for  om  <  a  <  o  m  , 
v(a)  -  0  for  o*m  <  a  <  , 

v(a)  >  0  for  <  a  <  , 


where  a  m  and  a  ^  are  the  unique  numbers  defined  by 


f3l(^*m)  “  m*  -  f31<ff*M)  ”  M* 


Clearly,  a  m  and  a  ^  satisfy 


am  <  a*m  ^  ac  ^  a*M  <  aH 


where  ac  is  the  Maxwell  stress  (Figure  1).  Clearly,  the  two  transformation 
stresses  EM  and  associated  with  shock  initiation  must  satisfy 


am  <  ^m  ^  a*m  ^  ac  ^  CT*M  ^  SM  <  CTM- 


(7.7) 


Property  (7.3)  of  V  is  readily  shown  to  imply  that  there  are  numbers  am' 
and  ,  with  am  <  am'  <  a*m  and  <  cr^,  and  such  that 


r(°)  <  cm(^  -  ^m)'1  for  am  <  °  <  <  \ 

r(o)  >  cM(oM  *  <7>'1  for  a  •  J 


(7.8) 


for  suitable  constants  cm  <  0  and  c^  >  0  .  A  schematic  graph  of  v(a)  is 
shown  in  Figure  7(b). 

If  <p  stands  for  the  inverse  of  the  restriction  of  v  to  (<rm , cr*m ]  + 
(<?  m,<7m)’  versi°n  (6.5)  of  the  kinetic  relation  may  be  put  in  the 
alternative  form 


o(t)  -  F(t)/A(s(t) )  -  <p(I(t)). 


Note  that  <p  has  a  discontinuity  at  the  origin  unless  a*m  -  ac  -  <7*^. 


(7.9) 


MiH 
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Suppose  for  the  moment  that  we  have  a  motion  taking  place  on  the  time 
interval  [tg,T]  and  involving  a  (3,1) -shock  located  at  x  -  s(t)  at  time  t. 
As  time  increases,  the  shock  location  s(t)  will  evolve  according  to  the 
kinetic  relation  (6.5),  which  in  present  notation  is 

s(t)  -  v(o(t))  -  v(F(t)/A(s(t))).  (7.10) 

Property  (7.8)  guarantees  that  the  moving  point  (s(t),F(t))  in  the 
(s,F) -plane  remains  in  the  set  corresponding  to  equilibrium  states  with 
(3 , 1) -shocks .  To  prove  this,  it  is  sufficient  to  show  that  the  stress  at 
the  shock  o"(t)  -  F(t)/A(s(t))  never  exceeds  and  is  never  less  than  am. 
Suppose  that,  at  the  initial  instant  tg,  one  has  o^'  <  o,(tg)  <  so  that 
the  second  of  the  inequalities  in  (7.8)  holds  at  time  tg.  We  shall  show 
that  <r(t)  <  Ofi  for  all  t  in  [tg,T].  Suppose  this  were  not  the  case.  Then 
there  would  be  instants  t  in  (tg,T]  at  which  cf(t)  -  aM;  let  t^*  be  the 
infimum  of  all  such  instants.  Clearly 

tg<ti*<T,  a(t)<aM  for  tg<t<t]*,  and  y(t]*)-oM,  (7.11) 

by  the  continuity  of  o(t) .  Now  let  tg*  be  the  supremum  of  the  set  of  all 
times  t  for  which  <F(t)  <  oM'  and  tg  <  t  <  t^*.  Then  for  tg*  <  t  <  t^*,  we 
have  Oft'  <  o(t)  <  so  that  (7.8)2  applies,  and  v(a(t))  >  0,  whence  by 
(7.10),  s(t)  >  0  during  this  time  interval  as  well.  It  follows  that  we  may 
express  t  as  a  function  of  s:  t  -  t(s),  and  thus  regard  <T(t(s))  -  <T(s)  as  a 
function  of  s  as  well.  Then  by  (4.17),  (7.10), 

o'(s)A(s)  +  <f(s)A' (s)  -  d/ds { F( t (s) ) }  -  F( t(s) )/S(t(s) ) 

-  F(£(s) )/v(o(s) ) .  (7.12) 

Let  X  -  max  | F( t ) | ,  tg  <  t  <  T,  be  the  maximum  loading  rate  during  the 
motion.  Then  by  (7.8)2,  (7-12)  yields 

o' (s) A(s)  +  o(s)A'(s)  <  (A/cM)(oM  -  o(s)),  sQ*  <  s  <  s^*,  (7.13) 

where  sQ*  -  s(t0*),  s^*  -  s(t^*).  Integration  of  this  linear  differential 
inequality  and  using  the  fact  that  A(sg*)  5  A(s)  leads  to 
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<f(s)  £  -  (<tm  -  <70)exp{ -(A/cm)Ja(x)  ‘^dx)  ,  sQ*  <  s  <  S3*,  (7.14) 

so* 

with  <7q  -  a(s0*).  In  particular,  this  gives  <7(s^*)  -  ^(tj*)  < 
contradicting  (7.11).  It  follows  that  “(t)  <  ojj  for  tQ  <  t  <  T.  Thus  the 
stress  <x(t)  at  the  shock  in  a  motion  governed  by  (7.10)  never  exceeds  a 
similar  argument  shows  that  “(t)  is  never  less  than  am. 

7 . 2  Macroscopic  response. 

We  now  elucidate  the  nature  of  the  macroscopic  response  of  the  bar 
under  various  force-controlled  programs  of  loading  for  the  class  of  kinetic 
relations  described  above.  First,  let  F(t)  -  At,  0  <  t  <  T,  corresponding 
to  loading  at  a  constant  rate  A  from  the  undeformed  state.  Assume  that  the 
final  value  of  force  is  such  that  F(T)  >  <*maM-  The  force-elongation 
response  answering  to  this  loading  is  shown  schematically  in  Figure  8(a). 
After  loading  begins,  the  point  (6(t),F(t))  rises  from  the  origin  0  along 
the  response  curve  0A0*  associated  with  smooth  fields  of  the  type  (4.1),  so 
that  S(t)  -  A^(F(t) ) ,  and  there  is  no  dissipation.  When  the  force  reaches 
the  level  corresponding  to  the  point  0*  in  Figure  8(a),  the  stress  at  x  -  0 
coincides  with  the  larger  transformation  stress  ZM,  and  a  (3,1) -shock  is 
created  at  x  -  0  according  to  the  shock  initiation  criterion  of  the  preced¬ 
ing  section.  The  kinetic  relation  (7.10)  takes  over,  and  the  initial  shock 
velocity  has  the  value  v(Zfl) .  By  (7.7),  5*  >  o*M,  so  that  by  (7.4),  v(Zj^)  > 
0.  If  Em  >  a  h*  then  v(Zfl)  >  0,  the  initial  shock  velocity  is  positive,  and 
there  is  a  discontinuity  in  slope  in  the  F-6  response  curve  at  0*,  as  shown 
in  the  figure.  If  Z^  -  the  initial  shock  velocity  is  zero,  and  the 

slope  of  the  F-6  curve  is  continuous  at  0*.  Under  the  control  of  (7.10), 

s(t)  increases  with  t,  the  elongation  is  given  by  5(t)  -  A31 (F( t) , s ( t) ) , 

and  the  point  (5(t),F(t))  moves  along  the  curve  0*BP,  both  F(t)  and  S(t) 
increasing.  This  stage  of  the  process  Is  accompanied  by  dissipation.  When 
the  force  has  increased  to  the  value  associated  with  point  P  in  Figure 
8(a),  the  shock  has  arrived  at  the  end  x  -  L  of  the  bar,  all  particles  are 

in  phase  3,  and  the  field  is  smooth.  As  the  force  continues  to  increase, 

the  response  is  that  of  smooth,  dissipation- free  phase-3  deformations,  with 
6( t)  -  A33(F(t)).  The  loading  terminates  at  time  T,  corresponding  to  the 
point  Z  in  the  figure. 

Now  suppose  that  F(t)  is  decreased  at  the  constant  rate  A  from  its 
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largest  value  F(T)  to  zero.  The  response  curve  at  first  follows  the  arc 
ZPQ*  corresponding  to  smooth  phase-3  fields,  and  the  response  is  dissipa¬ 
tion-free  ;  at  Q*,  the  stress  at  the  larger  end  x  -  L  of  the  bar  has  dimin¬ 
ished  to  the  smaller  transformation  stress  5^,  and  a  leftward  moving 
(3,1) -shock  emerges.  The  sign  of  v  in  (7.10)  is  now  negative,  s(t),  F(t) 
and  S( t)  all  decrease,  and  the  arc  Q*CA  is  traced  out  as  the  shock  returns 
to  x  -  0,  dissipating  as  it  moves.  As  F(t)  decreases  to  zero  from  its  value 
at  A,  S( t)  -  A^(F( t) )  again  along  the  arc  AO,  and  the  final  stage  of 
unloading  takes  place  without  dissipation.  The  total  dissipation  in  the 
loading  cycle  is  of  course  precisely  the  area  of  the  hysteresis  loop 
A0*PQ* . 

If  the  loading  rate  A  were  changed,  the  loading  and  unloading  "yield 
points"  0*  and  Q*  would  remain  the  same,  but  the  arcs  0*BP  and  Q*CA 
associated  with  the  dissipative  portion  of  the  process  would  change.  The 
macroscopic  response  is  thus  rate  -  dependent . 

Consider  now  a  modified  version  of  the  loading  history  described  above 
in  which  the  force  F(t) ,  after  arriving  at  the  value  associated  with  the 
point  B  in  Figure  8(a),  is  decreased,  and  then  ultimately  increased  again. 
Figure  8(b)  shows  the  resulting  macroscopic  response.  As  before,  the 
response  curve  is  the  arc  00*B  during  the  initial  loading  phase,  the  por¬ 
tion  0*B  being  dissipative.  When  F(t)  begins  to  decrease,  v(F(t)/A(s(t) ) 
remains  positive  at  first,  and  (7.10)  requires  s(t)  to  continue  to 
increase,  accompanied  by  dissipation.  During  this  stage,  &( t)  will  also 
continue  to  increase,  generating  the  arc  BC  of  the  response  curve.  At  the 
point  C,  the  stress  a  at  the  shock  has  dropped  to  the  value  a  so  that  by 
(7.4),  v(7)  -  0  at  the  corresponding  instant.  If  a  ^  >  a  m  in  (7.4),  and  if 
F(t)  continues  to  decrease  below  its  value  at  C,  <T(t)  will  remain  in  the 
range  for  which  v(<7(t))  vanishes,  so  that  s(t)  -  0  and  the  shock  remains 
stationary.  During  this  portion  of  the  unloading  process,  the  corresponding 
part  CD  of  the  response  curve  lies  along  a  curve  of  constant  s,  and  there 
is  no  dissipation.  If  now  the  force  F(t)  is  increased,  the  Initial  portion 
of  the  reloading  process  takes  place  along  DC  and  is  dissipation- free .  If 
the  force  ultimately  increases  sufficiently  to  raise  the  stress  at  the 
shock  to  a  value  greater  than  a  the  shock  resumes  its  motion,  dissipa¬ 
tion  will  begin  again,  and  the  response  curve  will  proceed  along  CE.  This 
particular  force  history  illustrates  the  occurrence  of  dissipation- f ree 
unloading  with  a  stationary  shock. 
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If,  during  unloading  in  the  last  program,  the  force  had  been  decreased 
sufficiently  below  its  value  at  D,  the  stress  at  the  shock  would  diminish 
below  the  value  o*m,  causing  v(ZT(t))  to  become  negative,  forcing  the  shock 
to  move  left.  Figure  8(c)  shows  the  macroscopic  response  curve  00*BFI  for 
such  a  force  history,  together  with  the  response  on  reloading.  The  arcs 
00*,  CE,  and  GH  correspond  to  dissipation- free  periods  during  the  quasi¬ 


static  motion. 


The  macroscopic  response  of  the  bar  during  the  loading  programs  just 
described  clearly  resembles  that  associated  with  visco-plastic  behavior  in 
several  respects.  One  feature  of  the  latter  kind  of  behavior  that  is  not 
present  here  is  that  of  permanent  strain.  By  abandoning  the  requirement  am 
>  0  in  (2.9)  and  thus  considering  a  stress-strain  curve  whose  local  minimum 
(Figure  1)  corresponds  to  a  compressive  stress  am,  one  can  introduce 
permanent  strain  into  the  macroscopic  response;  see  [18], 


The  form  of  the  kinetic  response  function  sketched  in  Figure  7(b)  sug¬ 
gests  consideration  of  the  limiting  case  in  which  the  function  < p  inverse  to 
v  is  a  step- function  as  shown  in  Figure  9: 


:D  -  {  * 

l  2. 


for  s  >  0, 
for  5  <  0, 


(7.15) 


where  Z^  and  are  the  shock  initiation  stresses,  am  <  5^  <  ac  <  Z^  < 
and  o c  is  the  Maxwell  stress.  One  shows  readily  that  the  macroscopic 
response  produced  by  the  kinetic  relation  is  rate  -  independent  and  is  of  the 
form  shown  in  Figure  9.  If  Zm  -  om  and  Z^  -  a the  quasi-static  motions 
permitted  by  the  kinetic  relation  are  maximally  dissipative  in  a  definite 
sense;  response  of  this  kind  is  discussed  in  detail  in  [17], 


Finally,  we  note  that  purely  conservative  (or  dissipation- free) 
response  of  the  kind  conventionally  associated  with  elastic  behavior 
results  from  choosing  the  inverse  kinetic  response  function  <p  to  be 


<p(  s)  -  <7C  for  -00  <  S  <  oo  , 


(7.16) 
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and  caking  both  shock  initiation  stresses  2^  and  Zq  to  be  equal  to  the 
Maxwell  stress  ac .  In  this  case,  the  macroscopic  response  is  independent  of 
past  history  and  is  as  shown  in  Figure  10. 
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FIGURE  I.  STRESS -STRAIN  CURVE. 


FIGURE  3.  EXAMPLES  OF  THE  FOUR  TYPES 


OF  TRANSITION  INSTANTS  t*.  IN  (4.7). 
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FIGURE  4.  ADMISSIBLE  DIRECTIONS 
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FIGURE  5.  SETS  <*  AND  ADMISSIBLE  DIRECTIONS  FOR  THE  UNIFORM  BAR 


FIGURE  7.  KINETIC  RESPONSE  FUNCTIONS  V  AND 
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cr  =  cr. 
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FIGURE  8(b).  FIRST  LOADING -UNLOADING- RELOADING  PATH 

ACCORDING  TO  KINETIC  RELATION  (7.10). 
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continuous  strains.  Boundary-value  problems  for  such  "unstable"  elastic  materials  often 
possess  an  infinity  of  solutions,  suggesting  that  the  theory  suffers  from  a  constitutive 
deficiency.  In  the  setting  of  the  one-dimensional  theory  of  bars  in  tension  ,  the  present 
paper  explores  the  consequences  of  supplementing  the  theory  with  further  constitutive  infor¬ 
mation.  This  additional  information  pertains  to  the  surface  of  strain  discontinuity  and 
consists  of  a  "kinetic  relation"  and  a  criterion  for  the  "initiation"  of  such  a  surface. 

We  show  that  the  quasi-static  response  of  the  bar  to  a  prescribed  force  history  is  then 
fully  determined.  In  particular,  we  observe  how  unstable  elastic  materials  can  be  used  to 
model  macroscopic  behavior  similar  to  that  associated  with  viscoolasticity . 
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